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Abstract 



I compute the hard spectator interaction amplitude in B — > 7r7r at NLO i.e. at 
,-Ch 1 O(cVg). This special part of the amplitude, whose LO starts at 0(a s ), is defined in 

the framework of QCD factorization. QCD factorization allows to separate the short- 
Q-ij and the long-distance physics in leading power in an expansion in Aqcd/ 77 ^; where 

(— i ' the short-distance physics can be calculated in a perturbative expansion in a s . 

In this calculation it is necessary to obtain an expansion of Feynman integrals in 
t^j- . powers of Aqcd/^- I will present a general method to obtain this expansion in a 

systematic way once the leading power is given as an input. This method is based on 
differential equation techniques and easy to implement in a computer algebra system. 



■ The numerical impact on amplitudes and branching ratios is considered. The NLO 

contributions of the hard spectator interactions are important but small enough for 
| perturbation theory to be valid. 

o 



1 Introduction 

In the last decades B physics has proven to be a promising field to determine parameters 
of the flavour sector with high precision. On the theoretical side QCD factorization [H [2] 
has turned out to be an appropriate tool to calculate B decay modes from first principles. 
Though the decay of the S-meson is caused by weak interactions, strong interactions 
play a dominant role. It is however not possible to handle the QCD effects completely 
perturbatively. This is due to the energy scales that are contained in the 5-meson: Whereas 
a s at the mass of the 6-quark is a small parameter, the bound state of the quarks leads to an 
energy scale of O(Aqqd) which spoils perturbation theory. The idea of QCD factorization 
is to separate these scales. At leading power in Kqcd/ttii, we obtain the amplitude for 
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Figure 1: Tree level, vertex correction and penguin contraction. These diagrams contribute 
to T l . 

B — > 7T7T in the following form: 



Two different types of quantities enter this formula. On the one hand the hadronic physics 
is contained in the form factor F B ^ n and the wave functions <Pbi and which will be 
defined more precisely in the next section. These quantities contain the information about 
the bound states of the mesons. They have to be determined by non-perturbative methods 
like QCD sum rules or lattice calculations. Alternatively, because they are at least partly 
process independent, they might be extracted in the future from experiment. On the other 
hand the hard scattering kernels T 1 and T n contain the physics of the hard scale 0(rrib) 
and the hard collinear scale O(yW^qcd) and can be calculated perturbatively. 

The Feynman diagrams that contribute to B — > titt can be distributed into two different 
classes. The class of diagrams where there is no gluon line connecting the spectator quark 
with the rest of the diagram (fig. [T]) contributes to T . We obtain T n by evaluating the 
hard spectator scattering diagrams, which are shown in LO in a s in fig. [21 The order ot 2 s 
corrections of T 11 are the topic of the present work. Through the soft momentum / of 
the constituent quark of the .B-meson the hard collinear scale */ 'Aqcd 77 ^ comes into play. 
This leads to the fact that in contrast to T 1 , which is completely governed by the scale ra&, 
T 11 comes with formally large logarithms. These logarithms cannot be resummed in the 
present QCD calculation. It will be shown by numerical analysis that the scale dependence 
of the hard spectator scattering amplitude and the absolute size of its NLO corrections are 
small enough for pertubation theory to be valid. 

My calculation of the hard spectator scattering amplitude is not the first one as it has 
been calculated recently by 0, H] • It is however the first pure QCD calculation, whereas 
[HI H] used the framework of soft-collinear effective theory (SCET) j5j El [7] an effective 
theory, where the expansion in Aqcd / n^b is performed at the level of the Lagrangian rather 
than of Feynman integrals. It is the main result of this paper to confirm the results of [Hill] 
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Figure 2: Hard spectator interactions at 0(a s ). This is the LO of T 



and to show by explicit calculation that pure QCD and SCET lead to the same result in 
this case. 

Hard spectator scattering corrections to the penguin diagram (third diagram of fig. [T]) 
are beyond the scope of this publication. This class of diagrams does not influence the 
cancellation of the scale dependence of the "tree amplitude" i.e. the diagrams of fig. [2] 
and higher order a s corrections. For phenomenological applications, however, the hard 
spectator penguin amplitudes, which have been recently calculated in [8], should be taken 
into account. Also the order o? s of T 1 is important for phenomenological applications. This 
calculation has been partly performed by [HI [TO] . There the complete imaginary part and 
a preliminary result of the real part of the amplitude is given. 

From a technical point of view the calculation in this paper consists of the evaluation of 
about 60 one-loop Feynman diagrams. The challenges of this task are due to the fact that 
these diagrams come with up to five external legs and three independent ratios of scales. 
In order to reduce the number of master integrals and to perform power expansions of 
the Feynman integrals, integration by parts methods and differential equation techniques 
will prove appropriate tools. They provide a general method to obtain higher powers of a 
Feynman integral once the leading power is given. 

The paper is organized as follows: I define my notations in section In section [3] I 
show how to get T 11 at LO. In section H] I present a method which uses differential equation 
techniques and allows for the extraction of higher powers of Feynman integrals once the 
leading power is given. Sections is dedicated to the technical details of the calculation. 
After some remarks how to evaluate the Feynman diagrams occuring at NLO I show how 
to deal with meson wave functions at NLO. Especially the correct treatment of evancescent 
structures will be explained in detail. In section [6] the analytic results of the calculation 
will be given, whereas section [7] will provide the numerical analysis. I end up with the 
conclusions. 
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2 Notation and basic formulas 



2.1 Kinematics 

For the process B — > irir we will assign the momenta p and q to the pions (fig. [2D which 
fulfil the condition 

p 2 ,q 2 = 0. (2) 

This is the leading power approximation in Aqcd/^6- Let us define two Lorentz vectors 
n+, n_ by: 

n% = (1,0, 0,1), r£ = (1,0, 0,-1). (3) 

In the rest frame of the decaying meson p can be defined to be in the direction of n + and 
q to be in the direction of n_. Light cone coordinates for the Lorentz vector z^ are defined 
by: 

-, z ± = (oyy,o) (4) 



" x/2 ' ~ V2 
So one can decompose z M into: 

p ■ q p ■ q 

such that 

z± ■ P = Z± ■ q = 0. (6) 

2.2 Colour factors 

In our calculations we will use the following three colour factors, which arise from the 
SU(3) algebra: 

Cn = \, C ^ = %^ and °g = N c , (7) 
where N c = 3 is the number of colours. 

2.3 Meson wave functions 

The pion light cone distribution amplitude 0^ is defined by 

(7r(p)|g(2) Q [...]g^O^|O) z2=o = ^(^75W^ 1 ^e^' 2 7r (x). (8) 
The ellipsis [. . .] stands for the Wilson line 

[z, 0] = P exp (J dt ig s z ■ A(zt) S j 
4 



(9) 



which makes (jHJ) gauge invariant. For the definition of the 5-meson wave function <f> B i we 
need the special kinematics of the process. Following [2] let us define 

*f(z,p B ) = (0\q,(z)[. . .]b a (0)\B(p B )) = J ^Le-^f^PB). (10) 

In the calculation of matrix elements we get terms like: 

^tr(„4(/)<M/)) = J-^J d*ze^tr(Am B (z)). (11) 

We have to consider only the case that the amplitude A depends on / only through / ■ p: 

A = A{l-p) (12) 
In this case we can use the 5-meson wave function on the light cone which is given by [2]: 



,.2 1—0 



m a {z)\..]b P mB{p B )) 



3 Hard spectator interactions at LO 

The effective weak Hamiltonian we deal with is given by [TT] : 



n cS = %v: d v ub [CxOx + c 2 o 2 ] + h.c, 



(13) 



'o 

where ^ ^ 

de*Bi(0 = l an d / d^ B2 (O = 0- (14) 
o Jo 

It is now straightforward to write down the momentum projector of the 5-meson: 
^lr{A{2l ■ p)Hl)) 
= ^tr^s + m B ) l5 J d£ ($ B1 (0 + rt + $B2(0)A(£m 2 B ) (15) 
At this point we give the following definitions 

? - f%m (is) 

X B Jo 4 



^ / f iB»e^(0- (17) 

m B Jo 4 



where 



01 = {dp)v_ A {pb) v _ A , 

2 = {d i p j ) v _ A {p j b i ) v _ A . (19) 

Explicit expressions for the short-distance coefficients Cj can be obtained from [11]. The 
decay amplitude of B — > tttt is given by 

A(B -> tttt) = (irirlHeslB). (20) 

For later convenience we define 

A(B -> tttt) = .A(£ -> tttt) 1 + .4.(5 -> tttt) 11 (21) 

where .A 1 (A 11 ) belongs to the first (second) term of (pQ). Because A 1 and .A 11 contain 
different hadronic quantities, the renormalisation scale dependence of both of them has 
to vanish separately. So we can set their scales to different values // and /i n . As in A 1 
there occurs only the mass scale mj we can set /i 1 = m&. In A 11 there occurs also the 
hard-collinear scale yA.qcD 77 ^- As we will see this scale is an appropriate choice for /i n . 

Because we only deal with the tree amplitude and do not consider penguin contractions 
only the matrix elements of the operators 0\ and O2 are taken into account. These oper- 
ators close under renormalisation such that the corresponding hard spectator amplitude is 
independent of the renormalisation scale. 

The decay amplitudes of B — > nn can be written in terms of Oj as follows [T2j : 



where 



and 



-A(B°^n + n-) = [A^ai + \' p (a\ + r£ag)] 
-V2A(B- -> tt-tt ) = A^(ai + a 2 )Anr 

= [-A> 2 + A;K + rX)]^ (22) 



2m 2 

r% ^ m b (fi)(rh u (fi) + rh d (fi))' ^ 



For the LO and NLO results of the a, I refer to |12j . We define analogously to (121j) 

Oi = + (24) 

where the labels T and 'II' refer to the contribution to A 1 and A 11 . 

The leading order of the hard spectator interactions which start at 0(a s ) is shown in 
fig. [2j The hard spectator scattering kernel T n , which does not depend on the wave func- 
tions, can be obtained by calculating the transition matrix element between free external 
quarks, to which we assign the momenta shown in fig. [2J The variables x, x = 1 — x, y, y = 
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1 — y are the arguments of T n , which arise from the projection on the pion wave function 
(JHJ). In the sense of power counting we count all components of I of C(Aq CD ), while the 
components of p and q are 0(rrib) or exactly zero. We define the following quantities 

(S LL, (25) 
p ■ q p ■ q 

We will see that in the end the dependence on 9 vanishes in leading power such that we 
can use (IT5|) . 

We consider the three cases B° — > 7r + 7r~, B° — > 7r 7r° and B~ — > 7r _ 7r°. In the case, that 
the external quarks come with the flavour content of B° — > tt + tt~, the LO hard spectator 
amplitude for the effective operator O2 reads: 

<L(£° - * + *~) = 

(d(xp)u(xp) u{yq)d{yq)\0 2 \d{l)b(p + q- /)) spcct . = 
4na s C F N c — l -^d{l)^d{xp) u{xp)f{l - j B )b(p + q-l) 

x £, m B 

d(yq) ( ^ - 4^7^ (1 " 76)«(i7?), (26) 

V y yy J 

where the quark antiquark states in the input and output channels of the matrix element 
form colour singlets. The subscript "spect." means that only diagrams with a hard specta- 
tor interaction are taken into account. The amplitude of 0\ vanishes to this order in a s . In 
the case of B° — > 7r°7r° we get the tree amplitude from the matrix element of 0\. The case 
B~ — > 7r _ 7T° does not need to be considered separately, because from isospin symmetry 
follows [T31CE2]: 

VlA(B- -> 7T-7T ) = A(B° -> 7T + 7T-) + A(B° -> ttV). (27) 

On the other hand the full amplitude is the convolution of T n with the wave functions, 
given by (JTJ). To extract T 11 from (l26i) we need the wave functions with the same external 
states we have used in ( 1261) . i.e. we have to calculate the matrix elements (EJ) and ( |T0i) . 
where the pion or 5-meson states are replaced by free external quark states. To the order 
O(a s ) we get 

4>®*f>W) s |cZ(z- ? )e-^'(«(^)cZ(M)l4^)<(0)|0)-^=o 

= 27rN c 6(y' - y)dp(yq)u a (yq) 
^Xp( x ') = 2ttN c S(x' - x)u p (xp)d a (xp) (28) 

= y'tfa + e iI '-* + (O|^(z) i 6 a (O) i |d(06(p + g-Z)),-, x= o 
= 2irN c 5(l'- -r)d p (l)b a (p + q-l) 
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By using 



A {1) 

^spect. 



dxdydl 4^l a A X )^-^'iv)^t\-y'( l ) TU(1 \ X ^yJ )a'a/3'P 1 'l 



(29) 



we finally obtain: 

T ( >(x, y, l~) 



Ana. 



1 



a' 7 



2/ 



-lulu (1 - 7s) 

yy 



(30) 



It should be noted that only the first summand of the above equation contributes after 
performing the Dirac trace in four dimensions. The second summand is evanescent. This 
will be important, when we will calculate the NLO corrections of the wave functions (see 
section IB~2l . 

If we plug the hadronic wave functions defined by (jHJ) and (1151) into (1291) i.e. we calculate 
the matrix element (1261) between meson states instead of free quark states, we get for the 
LO amplitude @: 



AD 

^spcct. 



ifJ B c F 



AN} 



Ana s / dxdyd^Bi(0 ( l ) n(x)(l) ir (y) 



1 



(31) 



Following (CQ) and the conventions of [3] we write our amplitude in the form: 

• i 



ecu = -im% I dxdyd^Tl I (x,y,C,)f B ^Bi(OMn(x)f n (j) n (y). 



(32) 



where in the case of B — > tt + ti we define 



A 
A 



spect.l 
spect.2 



and in the case B — > w w we define 



.4 
.4 



spect.l 
spect.2 



(C^spect. 

(Ox) spect. 



(Cl)spect. 
(^2) spect. 



(33) 



(34) 



Because we use the NDR-scheme which preserves Fierz transformations for 0\ and O2, Tj l 
has the same form for both decay channels. From (13"TT) and (13"2"1) we get: 



n(i) 



Ana, 



C f 



4iV c 2 £xym B 



T, 



ii(i) 



(35) 



Aspect 



is used for the matrix elements of the operators Oi between both free external quarks and 



hadronic meson states. It should become clear from the context what is actually meant. 
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4 Calculation of Feynman diagrams with differential 
equations 



In this section I will discuss the extraction of subleading powers of Feynman integrals with 
the method of differential equations [TU [THl HE]. This method will prove to be easy to 
implement in a computer algebra system. The idea to obtain the analytic expansion of 
Feynman integrals by tracing them back to differential equations has first been proposed 
in [T3j. This method, which is demonstrated in [T3j by the one-loop two-point integral 
and in [15] by the two-loop sunrise diagram, uses differential equations with respect to the 
small or large parameter, in which the integral has to be expanded. 

In contrast to [TH [To] I will discuss the case where setting the small parameter to zero 
gives rise to new divergences. In this case the initial condition is not given by the differential 
equation itself and also cannot be obtained by calculation of the simpler integral that is 
defined by setting the expansion parameter to zero. It is not possible to give a general proof, 
but it seems to be a rule, that one needs the leading power as a "boundary condition". 
An efficient way to calculate the leading power of Feynman integrals is provided by the 
method of regions [T71 HH1 El EHj , whereas the subleading powers can be obtained from a 
differential equation. In the present section I will discuss which conditions the differential 
equation has to fulfil in order for this to work. 

Although the examples I use below are taken from the present calculation, this method 
is very general and can be used in any case in which the expansion of Feynman integrals 
in small parameters is needed. 

4.1 Description of the method 

We start with a (scalar) integral of the form 

J( Pl ,...,p n ,,m 1 ,...,m n ) = | ^ )dDi \ Dn (36) 

where the propagators are of the form Di = (k + pi) 2 — m 2 . We assume that there is only 
one mass hierarchy, i.e. there are two masses m <C M such that all of the momenta and 
masses pi and are of 0(m) or of (D(M). We expand ( 1361) in ^ by replacing all small 
momenta and masses by Pi — > Xpi and expand in A. After the expansion the bookkeeping 
parameter A can be set to 1. 

We obtain a differential equation for / by differentiating the integrand in f)36p with 
respect to A. This gives rise to new Feynman integrals with propagators of the form 

i 

and scalar products k ■ p^ in the numerator. Those Feynman integrals, however, can be 
reduced to the original integral and to simpler integrals (i.e. integrals that contain less 
propagators in the denominator) by using integration by parts identities. 
Finally we obtain for fl36l) a differential equation of the form 

4rI(X) = h(X)I(X)+g(X) (37) 
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where h(X) contains only rational functions of A and g(X) can be expressed by Feynman 
integrals with a reduced number of propagators. It is easy to see that h and g are unique 
if and only if / and the integrals contained in g are master integrals with respect to IBP- 
identities, i.e. they cannot be reduced to simpler integrals by IBP-identities. If /(A) is 
divergent in e = ^p, /, h and g have to be expanded in e: 

/ = £v 

i 

h = J2 hiei 

i 

9 = (38) 

i 

Plugging (1381) into (l3"Tj) gives a system of differential equations for Ii, similar to ( 1371) . In 
the next paragraph we will consider an example for this case. 

First let us assume that h(X) and g(X) have the following asymptotic behaviour in A: 

h(X) = /i (0) + Xh w + . . . 

g(X) = ^TAy^(lnA) (39) 

3 

i.e. h starts at A , and we allow that g starts at a negative power of A. We count In A as 
C(A°) so the gW may depend on In A. This dependence, however, has to be such that 

limAc/ (j) (lnA) = 0. (40) 
The condition (j4T)l) is fulfilled, if the are of the form of a finite sum 

m 

^a n ln"A. (41) 

n=no 

The limit m — ^ oo however can spoil the expansion (1391) . E.g. e _lnA = -r so the condition 
( H0|) is not fulfilled, which is due to the fact that we must not change the order of the limits 
A — » and m — > oo. 

Further we assume that also /(A) starts at A 

/(A) = J (0) (In A) + A/ (1) (In A) + . . . (42) 

and plug this into (1371) such that we obtain an equation which gives recursively: 

A 4 J« = £ dX'X' 1 " 1 h^I^-V (In A') + (In A') J . (43) 

I want to stress that, because /i starts at 0(A°), ( )43l) is a recurrence relation, i.e. does 
not mix into 1^ if j > i. As the integral is only well defined if % > 1, we need the leading 
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power l' ) as "boundary condition" and (HHj) will give us all the higher powers in A. It is 
easy to implement (l4*H|) in a computer algebra system, because we just need the integration 
of polynomials and finite powers of logarithms. 
A modification is needed if h starts at A -1 i.e. 

h = -l + h (f>) + .„. (44) 
A 

By replacing / = X n I we obtain the differential equation 

*I=(" + h)l + \»g (45) 



d\ \\ 
which is similar to (l3?l) and leads to 

c?A'A' i+ri_1 ^ hWl^-fiQnX') + ^^^(ln A')J , (46) 

which is valid for i > 1 — n. So, if / starts at C(A~ n ), the subleading powers result from 
the leading power. 

4.2 Examples 

We start with a pedagogic example: 
Example 4.1. 

I=[^L- I (47) 

J {2n) d k 2 {k 2 - \)(k 2 - 1) y } 

where A< 1. The exact expression for this integral is given by: 

i In A i 



lnA(l + A + A 2 + ...). (48) 



(4vr) 2 1 - A (4tt) 2 

We see that / diverges for A — >• 0. As described e.g. in [20] we can obtain the leading power 
by expanding the integrand in the regions k ~ y/\ and k ~ 1 . This leads in the first region 
to 

rl d k —1 7 /1 \ 

" w;2 , ;2 rr =-7— ^-r(l + e) - + l-lnA (49) 
(27r) d F(P - A) (47r) 2 ~ e \ e / 

and in the second region to 

d d k 1 i . , (\ \ 

<^kW^T) = (^ T(1 + e) {-t + 1 ) (50) 

such that we finally obtain 

/(°)(lnA) = ^lnA. (51) 
11 



This is the result we obtain from the leading power of (I48p . We write the derivative of I 
with respect to A in the following form: 



dX 1 - A 



d d k 



(2n) d k 2 (k 2 — \y 



(52) 



We obtained the right hand side of (j52p by decomposing -^-1 into partial fractions. Of 
course this decomposition is not unique which is due to the fact that / itself is not a 
master integral but can be further simplified by partial fractioning. From (1521) and (1371) 
we get: 

h = — !— = 1 + A + A 2 + ... 
1 — A 

9 = (hw^ = wf (x ' 1+1+x+ --- ] (53) 



such that the coefficients in the expansion in A according to (I39|) do not depend on the 
power label (k): 

/i« = l and 9 {k) = j^y r (54) 
We obtain for the recurrence relation 

K» = l jf dYY" gj^OnVH^ . (55) 

Using the initial value (IBTj) it is easy to prove by induction 

I {k) (In A) = In A VA; > 0. (56) 

(An) 2 



This result coincides with ( 1481 ). 

The first nontrivial example, we want to consider, is the following three-point integral: 



Example 4.2. 

d d k 



(2n) d k 2 (k + un_ + l) 2 (k + n + + n_) 2 ' 



(57) 



Here n + and n_ are collinear Lorentz vectors, which fulfil = n 2 ^ = and n + -n_ = |, 
u is a real number between and 1 and I is a Lorentz vector with I 2 = and P 1. 
Furthermore we define 

i = 2l-n + and 6 = 2l-n- (58) 

We expand / in /, so we make the replacement I — > XI and differentiate / with respect to 
A. The integral is not divergent in e such that we obtain a differential equation of the form 
(13*7)) where the Taylor series of h(X) starts at A as in (1391) . In g(X) only two-point integrals 
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occur, which are easy to calculate. I do not want to give the explicit expressions for h and 
g because they are complicated, their exact form is not needed to understand this example 
and they can be handled by a computer algebra system. Because the leading power of / is 
of C(A°), (H31) gives all of the subleading powers. 

We obtain the leading power as follows: First we have to identify the regions, which 
contribute at leading power. If we decompose k into 

k" = 2k ■ n+nt + 2k ■ n^n% + k\_ (59) 

we note that the only regions, which remain at leading power, are the hard region k^ ~ 1 
and the hard-collinear region 

k ■ n + ~ 1 
k ■ n_ ~ A 

kl ~ V\. (60) 

The soft region k^ ~ A leads at leading power to a scaleless integral, which vanishes in 
dimensional regularisation. In the hard region we expand the integrand to 

(61) 



k 2 (k + unJ) 2 (k + n + + n_) 2 



By introducing a convenient Feynman parametrisation we obtain for the (4— 2e)-dimensional 
integral over (IBTj) : 

T(l + e) exp(iire)- ( Hl - U) - I ln 2 (l - u)) . (62) 



(4vr) 2 - e v ' ^ v 'u V e 2 
In the hard-collinear region we expand the integrand to 

1 



k\k + an„ + 6n + ) 2 (2k ■ n+ + 1) ' 

The integral over (16"3"1) gives: 



(63) 



r(l + e) exp(ivre)- ^ ^ + 2Li 2 (w) + - ln 2 (l - u) 



+ lnwln(l-u)+ln(l-it)ln^V (64) 

Adding f[6"2l and (16"4"|) together we get the leading power of (1371) : 

= -(2Li 2 (u) + lnwln(l - u) + ln(l - u)]n9) . (65) 
(47r)^ it 
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By plugging (l6"oT) into (HHI) we obtain I at 0(A): 
7« = 

1 

lnw 



i 1 
(47r) 2 it 



ln(l-u)ln0 ln(l-u)ln« , , 2Li 2 (wl 
2 + lnw + — - + — i - + ln£ + — 



it 



it 



a 



^ln(l — u) ln(l — u)\n6 ln(l — u) \nu ^ 



a 



1 — u u 
ln£ 2Li 2 (w 
1 — it u 

Now we want to consider the following four-point integral 
Example 4.3. 

d d k 1 



it 



(66) 



(2yr) rf fc 2 (fc + n_) 2 (A; + I - n + ) 2 (k + / - un + f 



(67) 



where we used the same variables, which were introduced in (1371) . This example is very 
special, because in this case our method will allow us to obtain not only the subleading 
but also the leading power in /. I is divergent in e such that we obtain after the expansion 
( 1381) a system of differential equations of the following form: 



d 

dX 



I-i 



HqI + g-i 



^-I = h I + h x I-x + g - 
dX 

It turns out that in our example h takes the simple form 

2 + 2e 



A 



(6* 



(69) 



such that analogously to 



we can transform 
d 



into 



d 



AV 



d\ 



(A 2 / 



-2A/_! + A V 



(70) 



This system of differential equations can easily be integrated to: 



1 -1 



rW 



1 



A* +2 Jo 
1 



d\'\' l+l g { l- l) 



\i+2 



d\'X' l+1 



-2/2 + gt 1] 



(71) 
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where the superscript (i) denotes the order in A as in (IHUj) and (f4"2"l) . Both I_ x and J start 
at (^(A" 1 ). Because ( I7T1) is valid for i > — 1, it gives us the leading power expression, which 
reads: 

= ?rW r ( 1 + e 4 C-- 1 -^-^) < 72 > 

(47r)^ e u£ \ e 1 — u J 

where £ = 21 ■ n + as in the example above. The exact expression for (ISTj) can be obtained 
from [2Tj . Thereby (1721) can be tested. 

In the last paragraph I want to return to Example 14.21 I will show how we can use 
differential equations to prove that the integral (157|) depends in leading power only on the 
soft kinematical variable 9 = 21 ■ n_ and not on £ = 21 ■ n + . We need derivatives of the 
integral with respect to £ and 9, which we have to express through derivatives with respect 
to P. These derivatives can be applied directly to the integrand, whose dependence on Z M 
is obvious. We start from the following equations: 







dl» 0£ Ol 2 ' 



which lead to 



"W 1 = kW-tfL + DlLl. (74) 



where we have set I 2 = in ( J74"l) . Using ( J74"l) we can show that in leading power (151 
depends only on 6* and not on £. So we can simplify the calculation of the leading power 
by making the replacement Z M — > 0n+. The proof goes as follows: From (1421) we see that 
the statement "J^ ^ does not depend on £" is equivalent to 

^I(£A,0A) = 0(A). (75) 

Using the first equation of (1741) we get 

£|:/(£A, 0A) = 0(A)I(£A, 0A) + 0(A). (76) 

Because we know (e.g. from power counting) that 7(£A, #A) starts at A , (1751) is proven. 
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Figure 3: Gluon self energy 




Figure 4: Abelian diagrams 
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Figure 5: Nonabelian diagrams 

5 Technical details of the NLO calculation 

5.1 Evaluation of the Feynman diagrams 

The diagrams that contribute to Tj at NLO are listed in fig. [3] - EJ The diagrams are 
evaluated in leading power in Aqcd/"V We do not use an effective theory like SCET but 
evaluate the diagrams in full QCD using free external quark states, to which we assign the 
momenta given in fig. [21 The Feynman integrals have to be evaluated in leading power in 
AQCT>/ m b- After reducing the number of Feynman integrals by integration by parts (IBP) 
identities (22l [23], we get the leading power of the master integrals using the method of 
regions (see e.g. [20J ) . In many cases the method of regions allows for a further reduction 
of master integrals. E.g. the first diagrams in the second line of fig. [5] comes with the scalar 
integral 



d d k 1 



(2vr) d k 2 (k + xp- l) 2 (k + xp + yq - l) 2 ((k + p + q - I) 2 - m 



(77) 



This integral can be calculated in leading power by setting 6 = 0, i.e. we make the replace- 
ment P — > £g M where £ and 9 are defined in (J25l) . This can be seen as follows: Counting 
soft momenta as O(X) and hard momenta as 0(rrib) the regions of space where (J771) gives 
a leading power contribution are 

fc" ~ A 

k ■ p ~ A k^_ ~ v^A k ■ q ~ nib- 

In these regions P occurs only in the combination I ■ p. So we can make the replacement 
Z p — > It can be easily seen that this replacement allows to reduce (1771) to three-point 
functions by decomposing the integrand into partial fractions. Alternatively one can use 
the exact expression (11 3 7ft for the four-point integral with one massive propagator line, 
which is given in appendix [A] After taking the leading power it can easily be seen that we 
get the same result as by just making the replacement Z M — > 

Tensor integrals of the form J -^^d(k^, k^k v ', ...) /Denominator can be reduced to scalar 
integrals using the methods of [21]. This reduction may make it necessary to calculate 
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subleading powers of Feynman integrals to obtain some diagrams in leading power. Fur- 
thermore, subleading powers of Feynman integrals are necessary in some diagrams, where 
the leading power vanishes because of the equations of motion. The methods of the last 
section, however, allow for an extraction of the subleading powers once the leading powers 
have been calculated. 

It is instructive and helps to avoid mistakes to obtain the Feynman integrals in two 
independent ways. Instead of using arguments depending on power counting we can obtain 
the power expansion of the integrals by calculating the exact integrals and expanding the 
results. The exact expressions of massless four-point integrals are given in [21]. Quite 
general expressions for four-point integrals with one massive propagator are given in ap- 
pendix |A] All of the integrals that were used for the present calculation have passed this 
independent test. 



5.2 Wavefunction contributions 
5.2.1 General remarks 

It has already been demonstrated in section [3] how in principle we can extract the scattering 
kernel T n of (CD) from the amplitude if we know the wave functions. T n does not depend 
on the hadronic physics and on the form of the wave function <f) n and <pB in particular, so 
we can get T 11 by calculating the matrix elements of the effective operators between free 
quark states carrying the momenta shown in fig. [2] on page [3J Because we calculate T n in 
NLO we need unlike as in section [3]the wave functions up to NLO. Let us write the second 
term of (CQ) in the following formal way: 

Apect. = <Pn ® fa ® <\>b ® T u . (78) 
All of the objects arising in (ITS]) have their perturbative series in a s , so ( 17H|) becomes 

# ® # ® <4 0) ® T n « (79) 
0« ® 0f ® <4 0) ® T n « + 0f ® 0« ® <pf ® T n « + 
® # ® 0? ® + # ® 0( Q) ® 0g } ® T"( 2 ) 



A {1) 

*^spcct. 

A {2) 

"^spcct. 



where the superscript (i) denotes the ordeiH in a s . In order to get T 11 ^ we have to calculate 
"4-spect.> </4 and (j>g for our final states. Then T 11 ^ is given by 

® 0f ® <4° } ® T 11 ^ = (80) 
4L. - ® ® 02 } ® - # ® 0« ® 0g) ® T 1 ^) - 
3 0g> ® T n « 



3 Note that the hard spectator scattering kernel starts at 0(a s ). So we call T n W the LO and T n ( 2 ) 
the NLO. 
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Figure 6: Example for diagrams which obviously belong to the form factor 



At this point a subtlety occurs. Let us have a closer look to the factorization formula (00). 
By calculating the first order in a s of the partonic form factor F B ^ n '^\ which is defined 
by free quark states instead of hadronic external states, we see that it can be written in 
the form 

F^)=0(°)®0g^7tLc,- (81) 

But T f ^ mfact is not part of T u . So we have to modify (180]) insofar as we have to subtract 
the right hand side of ( 1HTT) from the right hand side of (IBUl) : 



^ )®T II(2) = (82) 

# ® <j>™ ® <4 1} ® r 1 ^) - ® ® $ ® r^Ut. ® ^« 

In (1821) we did not include the term _F B_>7r >( 2 ) ® because it is obviously identical with 
the diagrams where the gluons do not interact with the emitted pion (e.g. those of fig. [6]). 
Those diagrams where not considered in the last section. So we do not have to consider 
them here. 

The wave functions for free external quark states are given at LO by (1281) . At NLO 
there exist three possible contractions: The two external quark states can be connected by 
a gluon propagator or one of the external quarks can be connected to the eikonal Wilson 
line of the wave function (fig. [7j). The diagrams of fig. [7(a), (b) and (c) give the order a s of 
the "pion wave function for free quarks", i.e. we have replaced the pion final state (7r(p)| in 
(IE]) by the free quark state (q'(xp)q(xp)\. The Fourier transformed wave function (x 1 ) 
is defined analogously to ( 1281) . For the diagrams in fig. [7J (a), (b) and (c) respectively we 
get: 



d d k S(x - x' - fi) - 5(x - x') 



-H + q'{xp) 



>SV) = ^ia,C F N c j ^ — 



q(xp)y 



xi 
1 



— X'i 



q' a (xp) 



a 
P 
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b( P - 



i-i) m 
(d) 



(e) 



(f) 



Figure 7: NLO contributions to the meson wave functions. The dashed line stands for the 
eikonal Wilson line which makes the wave functions gauge invariant. 



4 (c),(i). 



X 



8n 2 ia s C F N c 



d d k 5(x'-x+^] 
(2n) d k 2 



q{xp)^' 7 



xp — 



+ x 



-n T q'{xp) 



(83) 



5.2.2 Evanescent operators 

At NLO the convolution of the wave functions with the tree level kernel T 11 '^ 1 ' gives rise 
to new Dirac structures, which, however, can in four dimensions be reduced to the tree 
level Dirac structures. So we obtain the tree level Dirac structures plus further evanescent 
structures, which vanish for d = 4 but give finite contributions if they are multiplied by 
UV-poles. We define our renormalisation scheme such that we subtract the UV-poles and 
these finite parts of the evanescent structures. 

The tree level kernel ( 130]) contains two Dirac structures where the second one is evanes- 
cent (after the the projection on the wave functions). We write T II( ^ in the following 
form: 

T 1 ^ (*, y, 1-) - ^T'Wa - 75) ® ( ^ - (1 - 75) (84) 

xr v y yy J 

where the symbol (§) stands for the "wrong contraction" of the Dirac indices i.e. the Dirac 
indices are given by 

[rW ® r 3 ] Q ,_. = ryVJv, ( 85 ) 
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as in (I3(jp . The "right contraction" is defined by the symbol ® i.e. writing the Dirac indices 
explicitly 

[r 1 ® r 2 ® r3 ]« W7 = r^r 2 ^. (86) 

In d = 4 the wrong and the right contraction are related by Fierz transformations. It is 
convenient and commonly used to define the renormalised wave functions in terms of the 
right contraction, i.e. to define <^ en ' by renormalising the operator qiz^^^q' (0) instead of 
<j{z)pq'(0) a - This is why we define our renormalisation scheme such that only the UV-finite 
part of the right contraction operators remains: Using the notation of (I8"4"j) - (1861) we define 
the following operators: 

1 2i> 
O (x,y,r) = -^7—7^(1 -75)® 7m(1 -75)®— (1 -75) (87) 

Oi(x t y,r) = t^7 m ®7m(1-75)®— (1 -75) (88) 
t x y 

2 (x,y,r) = /r/®f(l-75)®^(l-75). (89) 
I x yy 

The matrix elements of these operators are defined analogously to ( 1291) : 

(Oi) = / dx'dy'dl'~ (j) naa '(x')(j)nl3/3'(y')(t>Byy'(l'~)Oi a/a /3'/3' r '' r (x',y', (90) 



Note that (0\ + 2 ) is just the convolution of the tree level kernel (|84[) with the wave 
functions. Furthermore by using Fierz identities it is easy to prove that we have in four 
dimensions 

(O ) = (O,) (91) 
(0 2 ) = 0. (92) 

So we define the following evanescent operators: 

Ei = 2 (93) 
E 2 = O 1 -O (94) 
If ~ (2- d) 2 

®2*f(l - 75) (95) 
£4 = — |p7V7 T ®7r7A7 !y (l-75)®^7M7v(l-75) (96) 

where we have defined E$ and E4 for later convenience. Using those operator definitions we 
define our renormalisation scheme such that we subtract the UV-pole of (Oq) and the finite 
parts of (Ei) i.e. terms of the form -^(E), where (E) is an arbitrary evanescent structure. 
It is important to note that we do not subtract IR-poles, because they depend not only on 
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the operator but also on the external states the operator is sandwiched in between. They 
have to vanish in (1821) such that the hard scattering kernel is finite. Finally we obtain the 
same result as if we had regularised the IR-divergences by small quark and gluon masses 
because the evanescent structures vanish in d = 4. The renormalisation scheme defined 
above is the same scheme that was used in [3]. 

In the next step we will calculate the convolution integral of T u '^ with the NLO 
wave functions given by (IHHj) . i.e. we have to calculate the renormalised matrix elements of 
Ox + C> 2 at NLO. 



5.2.3 Wave function of the emitted pion 

First we consider the renormalisation of the emitted pion wave function: Because the 
contribution of the wave functions (j^f 1 and (p£ does not change the Dirac structure of 
the operators, we do not need to consider evanescent operators when we calculate the 
diagrams of fig. 0(a), (b). So for the emitted pion wave function these diagrams give after 
renormalisation: 



//-Oren. , /rjren.x (l),(a),(b) _ ^ a s n 

{u ± + u 2 ; cmitted - -^t^f 



+ 21n^ * * (£>!) 
em Aim / y 

— + 2 ln/ii R ) (2 + In 3/ + In y) (0 2 > (O) 

^IR / 



(97) 



where the LO matrix elements (Oi)^ can be obtained from ( 1261) . Note that we kept the IR- 
pole times the evanescent matrix element (C^)^- This is needed for consistency because 
we also kept similar terms in the QCD-calculation of .Aspect.- Furthermore it allows us to 
show that all IR-divergences vanish. 

The diagram in fig. [7(c) mixes different Dirac structures. So we have to include evanes- 
cent operators in the renormalisation. In the case of the emitted pion wave function the 
operator 0\ does not mix under renormalisation with the evanescent operator E\ ( 193|) . We 
obtain for the renormalised matrix element: 

(OD^L = ~C F ^ (~ + 21n ^v) ( 0l >'»> (98) 

47r y v 6 ir a^ir / 

The matrix element of E\ however has an overlap with 0\. 



(^l)emittcd - — C F \Z — + 21n 



47T \ e uv em 



(-2y\ny-2y\ny) {E 1 ) 



(0) 



_ 4( |h» + ^!](o 1 )» 



(99) 

The renormalisation prescription tells us to subtract the UV-pole and the UV-finite part 
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of (Ei). So we obtain after renormalisation: 



(ED 



.\(l),(c) 
emitted 



-^ + 21n^X 



-2y\ny-2y\ny) (E^ 



(0) 



+ 4[lny+^)(0 1 )W 



(100) 



Note that the evanescent operator E\ leads to a finite term 4(lny + ^^ l )(Oi)^\ which we 
would have missed if we had just dropped the evanescent operators. 



5.2.4 Wave function of the recoiled pion 

In the next step we consider the NLO contribution of the recoiled pion wave function. As 
in the case of the emitted pion the diagrams fig. [7(a), (b) do not lead to a mixing between 
the operators. Therefore we get: 

ven.\(i).(»),(b) _ a s ~ 21nx + 4x 



(o\ cn - + o T 2 en ) 



recoiled 



47T 



X 



— + 21n^X 
em /^ir 



em 



(Oi) 



(0) 



21n/x m (0 2 ) 



(101) 



Other than in the case of the emitted pion the operators 0\ and O2 mix the spinors of 
the recoiled pion and the 5-meson. Therefore we have to work in the operator basis of Oq 
and the evanescent operators and define our renormalisation scheme such that the finite 
parts of the matrix elements of the evanescent operators vanish. The diagram fig. [7J(c) 
contributes to the matrix element of the renormalised operator O r cn ': 



(OD 



.\(l).(c) 
recoiled 



2^C F I - 
47r Vent 



21n- 



/iuv\ xlnx 



In the case of the evanescent operators we keep the IR-pole: 

1 a 



(Oo) 



(o) 



(ED 



(l),(c) 
recoiled 



2 4vr 



- + 2\n f , m )^(E 4 r 
em / x 



1 ETren.\(l)>(c) 



' recoiled 



xYax 



x 



<£ 3 > (0) 



(102) 

(103) 
(104) 



//'/-jrcn. 1 7-ircn. , rTiren. \ (l),(c) 
1 ' 2 /r 



At the end of the day we obtain a contribution from diagram fig. [TJc): 

v ~ I '' 2 fall* 

L , — 7^7 7 ®7 t 7a7 ® — : 

A*m/ art" V J 

KOi) (0) ' 



la s x In x 

Cp 

2 4?r x 



' recoiled 

— -21n^X 

e IR A*IR 



)7t7a7 



yy 



(105) 



The very complicated but also very explicit form, in which the above equation was given, 
is rather convenient, because the QCD amplitude -Aspect. ( see (1751) - (1501) ) comes with the 
same Dirac structure and cancels the IR-pole of (11051) . 
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5.2.5 Wave function of the B-meson 



The a s corrections of the wave function of the 5-meson are given by the second row of 
fig. fTJ For the diagrams (d), (e) and (f) respectively they read: 



'Bafl 



l7T 2 ia q N r C f 



lir 2 ia,N r C f 



d d k 5(l'- -I- - k-) - 5(l'~ - r) 
(2n) d Wk 1 

d d k 5{l'~ - I- - k~) - 5(1'- - l~) 



/ 



(106) 
b a (p + q-l) 



(27l) d 
1 



in 2 ia.N r C 



q(l)l 



1 



{+ /+ m b 
d d k 5(1'- -I 



k 2 k~ 
l~b(p + q - I) 



(2ir) d 



k 2 



/ 



m b 



ltib(p + q-l) 



In the case of the S-meson only the diagrams in fig. [7(d), (e) give rise to UV-poles. 
Those diagrams however do not lead to a mixing of Oq and the evanescent operators and 
we do not have to deal with evanescent operators. 

First let us have a look at the convolution integral which belongs to the diagram in 
fig- Elf) : 

(d + o 2 )^ m = 

(A-n)Ha 2 s N c C 2 F - [ d 
x J 



(107) 



x 



<?c7 



— 2 l%<irl U (l-ls) 



(2-n) d 2(k + l) -pk 2 

i~ /— ft + m b 



(k + l) 



k 2 - 2k ■ (p + q - I) 



2p~ 



iMe —Qn V -lf.lv 1(1 - l5)q' e 



yy 



where q c , qr and q^ are the spinors carrying the flavour quantum numbers of the light 
constituent quark of the S-meson, the recoiled pion and the emitted pion respectively. In 
leading power (11071) is identical to the contribution of the two diagrams shown in fig. El 
which is given by: 

d d k 1 



(4ir)HaiN c Ci 



(2n) d k 2 (k + I - xp) 2 



x qd T 



(k + l) 



^7^(1-75)^^ 



m b 



ytf+xji- ft- / y 

qe ~,<7 T T^lv - lv — 



k 2 -2k- (p + q-l) 
x$ — ft — / 



Irb 



(yq + xp — k — I) 2 



(yq + xp — k — I) 



;1» 



7s J q e 



108) 
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Figure 8: Two diagrams which correspond in leading power to the contribution of the 
5-meson wavefunction (I107p . 




Figure 9: a s contributions to the form factor 

In f llOSp the leading power comes from the region where k is soft. In this region of space 
the integrand gets the form of the integrand in (I107p . so both contributions cancel. As we 
did not include the diagrams of fig. M in the last section we can skip the contribution of 
(EEU7J) here. 

The remaining contributions are the diagrams in fig. [7(d) and (e). Together they read: 

(Oi + 2 )2 ),(d) ' (e) = (109) 
a 2 s N c C 2 F ^q c ^q' r q r Y(l ~ l 5 )bq e (—9^ ~ — lulu J q' e x 

£x \y yy J 

1 +21n^xV4 + 21nO-2f— + 21n^V4- — "21n 2 e 



5.2.6 Form factor contribution 

Finally we have to calculate the contribution of (1HT1) . It is given by 

Aormfact. = ^ <8> ^ ® ® T^act. ® = (110) 

^f (1) ' u UVQ)lu{l - l5 )q' e (yq) T« (y). 

The form factor f^ ,u is the a s correction of the matrix element 

{q' r {xp)q r (xp)\q rl v (l ~ ls)b\b(p + q - l)q c {l)), (111) 
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where T^(y) is given by 



T (1) (y) = -6 ( - + - 18 + 3 (— —\ny-in 



,9 21ny 
2Li 2 y -ln 2 y+— * 

2/ 



2/ 



(3 + 2ztt) lny - (y <-> y) 



We get Z^'" by evaluating the diagrams in fig. [9] and obtain finally: 

1 _ _ / / 
Aormfact. = o? a N cCp—q^ q% 7 At -7 I '(l - 75) - 7^(1 - 75)' 

g e7 ,(l-75)^T( 1 )(y). 



112) 



.r 



: 7, 6 



113) 



6 NLO results 

6.1 Analytical results for T^ 1 and T}} 

After the analysis of the last chapter we finally obtain the 0(a 2 ) results for the hard 
spectator scattering kernels Tl\ which are defined by fl3"2l . Those expressions appear in 
convolution integrals with wave functions, where x, y and £ are the integration variables 
as defined in ( l32l) . The ultraviolet divergences are renormalised in the MS-scheme. The 
infrared divergences drop out after subtracting the wave function contributions from the 
amplitude. The infrared finiteness together with the finiteness of the convolution integrals 
ensures that the framework of QCD-factorization works at this order in a s . 
The explicit 0(a 2 ) contributions for T^ 2 read (see next page): 
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The a 2 corrections of the hard spectator interactions have already been calculated in 
[Sill]. However both of these calculations have been performed in the framework of SCET, 
while my result is a pure QCD calculation. In order to compare (1 1 141) - ( TTT7T) to (3j H] we 
have to take into account the definition of Xb- The SCET calculation naturally uses the 
\b defined by the HQET field for the 6-meson, while I define by QCD-fields. Those 
two definitons differ at 0(a s ), which has been discussed in [25] . The difference in the 
logarithmic moments of the 5-meson wave function does not play a role, because these 
moments occur first at NLO. Using the results of [25] it is easy to figure out with the help 
of a computer algebra system, that flll4l) - flll7l) reproduce the results of [31 H]- 



6.2 Convolution integrals and factorizability 



By looking at the hard scattering kernels of (II 14H - ( TrTTj) it is not obvious that there remain 
no singularities in the convolution integrals over wave functions (1321) . It is however possible 
to perform the integration analytically, which proves the factorizabilty. 

Regarding the 5-meson wave function we will obtain the result in terms of the quantities 
Xb and X n , which are defined in ffToT) and (TTTT) . The 7r-meson wave function is given in 
terms of Gegenbauer polynomials: 



\x 



Qxx 



n=l 



'1181 



Due to the symmetry properties of the pion the first non vanishing moment is a\. We 



30 



neglect for n > 2 and using f[3"2l we get for the NLO of A spcct : 
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The finiteness of the above equations proves factorization of the hard spectator interactions 
at NLO. 

Including the contributions of A s ' t l and A s ' t 2 the quantities aiji and a 2 ,n defined 
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CKM-parameters 

Kd m Kd Kb m I Kb/Kb I [21 7 

0.974 -0.23 0.041 0.09±0.025 (70±20)de£ 

Parameters of the 5-meson 

m B f B [27] jM- m Ai [3] A 2 [3] r B± r B o 



/ 

5.28GeV (210±19)MeV 1.56 ±0.17 -3.2 ± 1 11 ±4 1.67 ps 1.54 ps 

Parameters of the 7r-meson 

/f ff [291 EDI EI] U m w a\ a\ [321 EH] 

0.28 ±0.05 131MeV 130MeV 0.3 ± 0.15 



Quark and W-boson masses 
m b (m b ) m c (m b ) m t (m t ) [12] M w 

4.2 GeV (1.3±0.2) GeV 167 GeV 80.4 GeV 

Coupling constants 
225 MeV 1.16639 x 10~ 5 GeV" 
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Table 1: Input parameters, which were used in the numerical analysis. All parameters given 
without explicit citation can be found in [34]. Unless otherwise stated scale dependent 
quantities are given at fi = lGeV. 
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7 Numerical analysis 
7.1 Input parameters 

For my numerical analysis I use the parameters given in table [H The decay constant 
and the ratio jb^ Xb have been obtained by QCD sum rules in [27] and [28] respectively. 

The logarithmic moments Ai and A 2 where calculated in [3] using model light-cone wave 
functions for the B-meson [331 ESI E7J EH]- For the form factor f^ 71 I use the value from 



32 



[3"T] . which has been obtained by QCD sum rules. This value is consistent with quenched 
and recent unquenched lattice calculations [2HJ ED] • The first Gegenbauer moment of the 
pion wave function is zero due to G-parity while the second moment has been obtained by 
lattice simulations [321 133] . 

7.2 Amplitudes a\ and 

The QCD amplitudes a\ and a 2 are defined in [12]. Their hard scattering parts a^n and 
a2,n, i.e. the parts of a\ and a 2 , which contribute to A 11 (see lj2"Tl)). are plotted in fig. [TO] 
as functions of the renormalisation scale \x. The strong dependence on fi of the real part 
of LO is reduced at NLO. Taking the twist-3 contributions into account does not increase 
the /i-dependence too much. The imaginary part, which occurs first at NLO, is strongly 
dependent on the renormalisation scale. An appropriate choice for the scale of the hard 
scattering amplitude is the hard collinear scale 

/i hc = 1.5GeV (123) 

In the following numerical calculations we will evaluate a^u and 02,11 a t /-*hc- The vertex 
corrections A 1 will be evaluated at 

fi b = 4.8GeV. (124) 
Using the parameters of table [1] we obtain 

a x =1.015 + [0.039 + 0.018i] y + [-0.012] tw3 + [-0.029] LO 

+ [-0.010 - 0.03H] NLO 
a 2 =0.184 + [-0.171 - 0.080z] y + [0.038] tw3 + [0.096] LO 

+ [0.021 + 0.045z] NLO . (125) 

These equations are given in a form similar to (61) and (62) in [3]. The first number 
gives the tree contribution, the vertex corrections are indicated by the label V, the twist- 
3 contributions are labelled by tw3. The hard scattering part is separated into LO and 
NLO. The hadronic input parameters I used are slightly different from [3] and in contrast 
to [3] I evaluated all quantities, which belong to the hard scattering amplitude, at the hard 
collinear scale /ihc- This is why the values I get for a\ and a 2 are different from [3]. 

The hard scattering amplitudes a± u and a 2 ji together with their numerical errors read: 

a i n = - 0.051 ± 0.01 l(param.)+°;^6 ( sca ie) ± 0.012(tw3) 

+ [-0.031 ± 0.008(param.)l£°g (scale) ± 0.012(tw3)]z 

a 2 ,n =0.15 ± 0.03(param.)+°;^(scale) ± 0.04(tw3) 

+ [0.045 ± 0.012(param.)+°;^°(scale) ± 0.038(tw3)]z. (126) 

The first error comes from the error of the input parameters in table [U The scale uncer- 
tainty is obtained by varying fi^ between lGeV and 6GeV. The error labelled by tw3 gives 
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Figure 10: Contribution of the hard spectator corrections to a\ and 02 as a function of the 
renormalisation scale //. The upper two figures show the real part, where the LO is given 
by the dashed line, while the sum of LO and NLO is shown by the thick solid line. The 
twist-3 corrections are included in the graph given by the thin solid line. The third figure 
shows the imaginary part, which occurs first at 0(a 2 s ). So no distinction between LO and 
NLO is made. 
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the error of the twist-3 contribution. Within the scale uncertainty (I126p is compatible 
with The result I obtained in QCD comes with formally large logarithms lnAqcD/ m b- 
Without resummation these logarithms might spoil perturbation theory. However the er- 
ror coming from the scale uncertainty in (1126j) as well as the relative size of the NLO 
contributions are small enough for perturbation theory to be valid. 



7.3 Branching ratios 

The dependence of the CP-averaged branching ratios on the hard collinear scale is shown 
in fig. HH It is obvious that the NLO corrections reduce this dependence significantly. 
From the parameter set table [T] we obtain the following CP-averaged branching ratios 

10 6 BR(5 + -> vr+vr ) = 6.05l^(had.)^;|°(CKM)l^(scale) ± 0.27(sublead.) 
10 6 BR(B° -> tt+tT) = 9.41+|^(had.)^;°°(CKM)^(scale)+J;^(sublead.) 
10 6 BR(5° -> 7r°7r°) = 0.39^;^(had.)l°;?°(CKM)+°;^(scale)+°;^(sublead.). (127) 

The origin of the errors are the uncertainties of the hadronic parameters and the CKM 
parameters, the scale dependence and the subleading power contributions, i.e. twist-3 and 
annihilation contributions. The error arising from the scale dependence was estimated by 
varying fi b between 2GeV and 8GeV and fi hc between lGeV and 6GeV. If we compare 
(11271) to the experimental values 



10 6 BR(5 + -> 7r + 7T°) = 5.5 ±0.6 
10 6 BR(£° -> tt+O = 5.0 ±0.4 

10 6 BR(5° -»• ttV) = 1.45 ±0.29 (128) 

we note that BR(5 + — > tt + tt°) is in good agreement with the data. For B + — > n + ii 
and B° — > 7r + 7r _ QCD- factorization is expected to work well, because at tree level Wilson 
coefficients occur in the so called colour allowed combination C\ + C2/ N c ~ 1, while B° — ► 
7r°7r comes at tree level with C2 + C\/N c ~ 0.2 such that subleading power corrections 
are expected to be more important. On the other hand there are big uncertainties in 
the parameters occurring in the combinations iK^I/jj 57 ", jmrj^ and aj. In |3Dj and [3j 

these parameters were fitted by the experimental values ( 1128j) of BK(B + — > 7r + 7r°) and 
BR(B° -> 7T+7T-). Setting 

aJ(lGeV) = 0.39 (129) 



leads to 



\V ub \f»* ^ 0.80(1^1/+") 



default 



2.89 ( -p— ) . (130) 



U A B \J+ A B/ default 
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Figure 11: CP-averaged branching ratios as functions of the hard collinear scale /Xhc in units 
of 10~ 6 . In the graph with the dashed line only the leading order of the hard spectator 
scattering is contained, while in the solid line hard spectator scattering is taken into account 
up to NLO. 
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This leads to the following branching ratios: 

10 6 BR(£ + -> tt+tt ) = 5.5 ± 0.2(param.)i°j(scale) ± 0.6(sublead.) 
10 6 BR(5° -> tt+tT) = 5.0lg;|(param.)18;|(scale)lg;g(sublead.) 
10 6 BR(fl° -»■ ttV) = 0.77 ±0.3(param.)io;3 (scale) igj(sublead.). (131) 

The uncertainties of the quantities that occurred in (I129p and (11301) have not been con- 
sidered in the estimation of the errors in f l 1 3 1 f) . The B° — > 7t°it branching ratio obtained 
in (11311) is compatible with the value obtained in [3J. Though it is too low, due to the 
theoretical and experimental errors it is compatible with (11281) . 

8 Conclusions 

QCD factorization has turned out to be an appropriate tool to calculate B decay modes 
from first principles, because it allows us to disentangle systematically the perturbative 
physics and the non-perturbative physics. The present calculation showed that the hard 
spectator scattering amplitude factorizes up to 0(a 2 ), i.e. all infrared divergences cancel 
and there are no remaining endpoint singularities. The former point is obvious after the 
explicit calculation of T 11 and the latter point was shown by evaluating the convolution 
integral ([1]) analytically. The explicit expressions for the hard spectator scattering kernel 
(jll4j) - (H17p confirmed the result of [3J H]. So they are also a confirmation that the leading 
power of the amplitudes can be obtained by performing the power expansion at the level 
of Feynman integrals rather than at the level of the QCD Lagrangian using an effective 
theory like SCET, which was done in [31 Hj. 

The main challenges in the evaluation of Feynman integrals were due to the fact that 
the Feynman integrals came with up to five external legs and three independent rations of 
scales. The calculation of the Feynman integrals was made possible with the help of tools 
like integration by parts identities and differential equation techniques: In section H] it was 
shown how to get the expansion of Feynman integrals in powers of AqcdA^ by differential 
equations once the leading power is given. 

Because next to the m;,-scale also the hard-collinear scale ^/Aqcd^& enters the hard 
spectator scattering amplitude, large logarithms could spoil perturbation theory. However 
the numerical analysis showed a strong reduction of the scale uncertainty of T 11 . It also 
confirmed the observation of [3J, that the NLO of T 11 is numerically important but small 
enough for perturbation theory to be valid. 

Finally it is important to note that in the present calculation the contributions of 
penguin contractions and the effective penguin operators were not considered. Actually 
they play a dominant role in the branching ratios of B — ► Kir and CP asymmetries of 
B — > -Kir and should be taken into account in phenomenological applications. They have 
recently been published in [8]. Also the 0{a 2 s ) corrections of T 1 were not part of the present 
work. These contributions have been calculated in [§1 HD] • 
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Figure 12: Basic one-loop four-point intergral. The massive line, which carries the mass 
m, is indicated by the thick line. 
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A Massive four-point integral 

We consider the following massive four-point integral in d = 4 — 2e dimensions (fig. [T2l : 




(132) 



where 



D 2 
D 3 
Da 



D 



k 2 + ii] 

(k + pi) 2 + irj 

[k + pi +p 2 ) 2 + if) 

(k + pi+p 2 +P3+ Pa) 2 -m 2 + 



(133) 



Following [21] we introduce the 



external masses 



p 2 = m 2 (2 = 1,2,3,4) 



(134) 



and the Mandelstam variables 



s = (pi+P2) 2 , t = (p 2 +p 3 ) 2 - 



(135) 
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Furthermore we consider only the case, where 



m\ = and ml = m 2 . 



(136) 



The integral (11321) can be evaluated using the method of [2"T] . This paper gives explicit 
expressions for massless one-loop box integrals. It is however possible to extend the single 
steps of this paper to our case. 

So finally we obtain: 



2\<E 



hipi,P2,P3,PA) = h(s,t,m 2 1 ,m 2 3 ,m 2 '' 
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where / 



s+t—rnf—m^ 



m 2 (rrij — s)+st— rn'frn'^ ' 

The case m 2 — gives rise to further divergences and has to be considered separately: 
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where f m = gL^I. 
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